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O^l ■ We have observed the S- and Q-branch Raman spectra of H2, HD, and D2 adsorbed at 85 K and 

5.^ ' pressures up to 8 atm on single- walled, carbon nanotubes (SWNT). Comparative data for H2 on 

Oh. graphite and C60 were also collected. For each adsorbate, we observed a small shift in the Q-branch 

1 frequencies relative to the gas-phase values. To aid in interpreting this result, we constructed an 

H2-surface potential, including van der Waals and electrostatic terms. Computed shifts based on 
this potential are in good agreement with our data. 
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Single-walled, carbon nanotubes (SWNT) are nanoporous, all-surface macromolecules which may be ideal storage 
media for H2 [0J|]. In addition, H2 has been predicted to exhibit novel, quasi-lD phases in SWNT ropes H§J. 
Desorption experiments suggest the presence of a high-energy binding site for H2 in purified and sonicated SWNT 
material, which has been proposed to be indicative of charge transfer |5|||. On the other hand, transport data for 
SWNT suggest that H2 physisorbs |7]] . An improved understanding of the adsorption of H2 on SWNT is the motivation 
for the work described here. We have collected the S- and Q-branch Raman spectra of H2, HD, and D2 adsorbed 
[ on SWNT at 85 K and at pressures up to 8 atm. The magnitude of the shifts in the Q-branch frequencies of these 
species when adsorbed on SWNT is a measure of the strength of the adsorption potential. For physisorption, the 
shifts should be small, while in the event of charge transfer, the shifts should be large: the change of the H2 stretching 
frequency with the molecular charge is ^2000 cm -1 per electron |fj"5|| . The approach taken in this work is similar to 
that described in Ref. || , except that we have employed Raman rather than IR spectroscopy. 
t-H , Using a two-step, oxidative and acid reflux technique, we obtained SWNT samples containing minimal amorphous 
and multishell carbon (<15 wt-%) and residual metal catalyst (<10 wt-%, <2 atom-%), as determined by temperature 
programmed oxidation. This material was annealed in a vacuum of better than 1(T 7 Torr at 1200 K for 24 hours; such 
high-temperature "degassing" was recently shown to remove the carboxylic and related functional groups occluding 
the SWNT ends and side walls . The resulting material was pressed against indium foil and inserted in an optical 
pressure cell, which was placed in a cryostat (Janis Corp. Model VPF 700, 77-700 K) with a gas line to the cell. The 
—I ' cell was valved to a diffusion pump and to the cylinders of H2 (ultra-high purity, 99.999%, MG Industries), HD (97.4 
mol-%, with 1.4% H 2 and 1.2% D 2 , Isotec Inc.), and D 2 (ultra-high purity, 99.999%, MG Industries). In separate 
experiments, samples of Ceo powder (Alfa Aesar, 99.9%, packed under Ar) and freshly cleaved, pyrolytic graphite 
(Alfa Aesar) were loaded into our cell. 

Raman spectra were collected in the backscattering geometry, with A a i r =514.53 nm excitation provided by a mixed- 
gas, Ar/Kr-ion laser (Coherent Inc., Innova Spectrum). The incident beam was focused through the windows of the 
cryostat and pressure cell onto the sample surface. Backscattered Rayleigh light was rejected by a holographic "super- 
notch" filter (Kaiser, Inc.), and Raman light was focused through a 25 fim slit into a single-grating monochromator 
q . (Instruments S.A. HR460, grating: 1800 grooves/mm) equipped with a CCD. To correct our spectra for slight, in- 
strumental nonlinearity and to calibrate the spectra, one Hg line and several pairs of Ne lines close to the Q-branch 
lines of H2, D2, and HD were used. The Ar, Hg and Ne frequencies were taken from Refs. |l(| and |l7| . 

For all data reported here, the sample temperature was 85±5 K, and the pressure of the adsorbate gas on the 
sample was increased from vacuum to as high as 8 atm, in the case of H2 and D2, and to 5 atm in the case of HD. 
Typical S- and Q-branch spectra at 8 atm are shown in Figs.l and 2; we use the notation S(J) or Q(J) to denote 
transitions, in which J is the initial state. The frequencies, v a ds, of the observed S- and Q-branch transitions are given 
in Table 1; these values were determined by fitting Voigt-profile functions to the spectra using Humlfcek's algorithm 
Also provided in Table-1 are the frequencies, ^/ ree , of these transitions for free molecules, from Ref. Jig] . 
Our S-branch data suggest freely rotating molecules in the gas phase, and are in good agreement with recent inelastic 
neutron scattering data Eg]. Small, downshifted shoulders arc evident in the S-branch lines, but our experimental 
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resolution is insufficient to quantify these effects. Our Q-branch data for HD and D2 contain prominent bands which 
are slightly upshifted from the free-molecule frequencies (see Table 1), but, due to the small splitting of the J=0 
and J=l bands, the peak substructure could not be resolved. Our data for H2 in SWNT are considerably more 
illuminating, and are shown in expanded form in Fig. 3 at two pressures, along with comparative data for Cgo and 
highly-oriented pyrolytic graphite (HOPG) collected at 85±5 K and 4 atm. The best fits (solid curves) to the data 
in Fig. 3 were computed using a sum of Voigt components (dashed curves). 

The surface area of HOPG is extremely small, so the signal strength from adsorbed molecules should be negligible. 
Therefore, we assign the bands observed at 4161.3 and 4155.4 cm" 1 to gas-phase H2 near the surface; these frequencies 
are in reasonable agreement with the those reported by May, et al Jl3| |. The gas-phase components should also be 
present in the data on Cgo and SWNT, and so our analysis of these data proceeded by inserting components of 
frequency and width equal to those observed on HOPG, and then introducing separate Voigt components to describe 
the additional bands. Consistent with this approach, we noted an approximately twofold increase in the relative 
intensities of the bands attributed to gas-phase species, as the pressure was increased from 4 to 8 atm. In our analysis 
of the data on SWNT at 8 atm, the positions of the bands attributed to the gas phase were allowed to vary, and the 
frequencies were found to be ~0.1 cm -1 lower, again consistent with the data of May, et al fill . 

In addition to the gas-phase components, it is clear that the data for Cgo and SWNT contain additional bands 
which are lacking in the data on HOPG. We ascribe these components to adsorption in the interstices or internal 
pores in the respective lattices. Neutron powder diffraction work by FitzGerald, et al., has indicated the existence 
of an adsorbed phase in the octahedral sites of /cc-Cgo El]]; presumably, these sites are partially filled under the 
conditions of our study and account for the lower-frequency component present in our Q-branch data for H2 on Cgo- 

A reduction in the number of fitting parameters was achieved by setting equal the linewidth parameters of equally- 
spaced pairs of Q(0) and Q(l) peaks. All parameters were then refined until \ 2 was minimized and the frequencies 
were obtained. Our confidence in the number of Voigt components introduced for each data set is affirmed by the 
relative intensities between "pairs" of Voigt components: Iq(i)/Iq(o) =2.2-2.5, in good agreement with the equilibrium 
ortho-para ratio at 85 K. Also, the Q(1)-Q(0) spacing we find between these pairs (5.8-6.0 cm" 1 ) is consistent with 
the value (~5.97 cm" 1 ) reported by May, et al, under similar conditions 

Based on shifts observed in the tangential Raman bands of SWNT material previously dosed with H2, partial 
electron transfer from SWNT to the H 2 adsorbate has been proposed p|. Our experiment is uniquely suited to 
looking for such effects, given the extreme sensitivity of the H2 stretching frequency on the molecular charge state. 
Contrary to reports on certain zeolites and oxides S, no strongly shifted Q-branch lines were observed on SWNT. 
This would seem to rule out significant charge transfer under our experimental conditions. 

To aid in interpreting the small shifts found in our Q-branch data, we have constructed an interaction potential 
for H2 with graphene (a single sheet of graphite) and estimated the frequency shifts in two types of adsorption sites. 
For the case of H2 adsorbed on graphene, the potential is a sum of C-H van der Waals interactions, Ulj, and the 
electrostatic interaction, U e i, of the H2 static multipole moments with the static screening charges induced on the 
graphene surface Jl9| . Because the low- frequency dielectric response of graphene is metallic [go) , the screening charges 
may be represented by full image charges displaced from the graphene sheet by a distance zq. The electrostatic 
interaction can be expressed in terms of the H2 quadrupole (O) and hexadecapole ($) moments ]l9| ]: 
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where z c is the perpendicular distance from the molecular center to the surface atoms, zi = z c — zq is the location 
of the image with respect to the graphene surface, 9 is the angle of the molecular axis with respect to the surface 
normal, and k is the usual electrostatic constant. In view of the small surface corrugation, the holding potential is 
V(z c ,cos9) = U e i + Ulj(zi) + Ulj(z%), where the zi are the perpendicular distances of the H atoms relative to the 
surface, and the Ulj(zi) are the corresponding, pairwise van der Waals interactions. The C — H interactions are 
assumed to be of Lennard-Jones (12-6) form, and for simplicity, we smear out the C atoms along the surface. In this 
approximation, 

EM*) = Meo-Boh-H [^(^r ~ C-^f] , (2) 

-0 2- Z - 

where i? = 0.38 A -2 is the surface density of the C atoms. We adopted the values ec-H — 2.26 meV , oc-h — 2-76 A 

The least well-known parameter in our potential is zq [ p2[ . To obtain this parameter, we computed the minimum 
of the C-H 2 potential at the equilibrium distance z eq , where H 2 is preferentially oriented flat against the surface. 
In the case of a graphene sheet, the isotropic C-H2 holding potential has the same functional form as Eqn.(||), with 
a minimum at z eq = ac~H 2 - Using eo_H 2 =42.8 K and ac-H 2 — 2.97A J24|], the resulting well depth is 46.6 meV. 
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The values for the H2 moments, taken from Ref. p5[ , are — 0.48917 a.u.; $ = 0.230063 a.u.. These values yield 
zq = 1.36 A, a value consistent with previous ab initio and empirical estimates for this system |23]| . 

For the study of H2 in the ICs, we adapted a potential, U^j, already developed for a single nanotube p6| . The 
dispersion part of this potential ignores many-body effects, considers the C atoms to be smeared on the surface, and 
is isotropic. When averaged over the azimuthal and longitudinal coordinates, this potential takes the form p6| : 

Uij(z) = 3^ec-HC7 2 C - H [^(^) 10 hi(x)M n (x) (^L) A h {x)M,{x)] , (3) 

in which z is the distance from the axis of a nanotube with nuclear radius R; we define x — R/z and fi(x) = (R/z) 1 
, with the I being positive integers. The Mi are integrals defined in Ref. |^6|. The potential for an IC is obtained by 
summing over an assembly of three nanotubes and azimuthally averaging the result. 

The electrostatic portion of the potential in the ICs contains two contributions, denoted U\ and U2, which arise, 
respectively, from the interaction of the H2 quadrupole moment with the local electrostatic field of the SWNT, and 
the interaction of the H2 static multipole moments with the screening charges induced on the surface. U\ is given by 
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in which the 0jj are the components of the quadrupole moment and <E> is the local electrostatic field, which was 
calculated from first principles for "zigzag " SWNT with radius R = 6.9 A (27). To evaluate U 2 , we included the IT 
quadrupole moment (which is the first nonzero permanent moment of an axial quadrupole), and the hexadecapole 
moment. This problem is quite complicated, so we introduce three simplifications. Due to the small size of the 
H2 molecule, we consider its multipoles to be "point"-like. Next, to simplify the complex geometry of the IC, we 
represent the SWNT walls with graphene planes. The locations of the image charges are determined by the boundary 
condition requiring that the potential be zero at each plane. Finally, we assumed that the graphene planes are perfectly 
conducting, so that the static and dynamic responses are represented by images. We adopted the same value for zq 
and the same form for the interaction as in the previous case. 

To compute the vibrational frequency shift, the total holding potential U must be expressed as a function of the 
relative bondlength change, £ = (r — r e )/r e , relative to the equilibrium value, r e , and also of the configuration r of 
the molecule. Given a sufficiently small H2-substrate interaction, one can treat U and the anharmonic terms in the 
potential energy function of the free H2 molecule as perturbations to the harmonic oscillator Hamiltonian ||. Using 
first- and second-order perturbation theory, the change in the frequency of the fundamental (n=0^n=l) transition 
due to U is [||||§: 

B B 
Auj = Auj^o = — - < U" - 3aU' > T +0(—^) 2 . (5) 

nuj e ncj e 

in which B e is the equilibrium rotational constant, a is the anharmonicity, primes denote derivatives with respect to 
£, and the angled brackets denote the average over all r. The expressions for 0(r) and 3>(r) were derived via a model 
in which the electron cloud is concentrated between the nuclei rather than about each individual nucleus. Assuming 
an H2 bondlength of 0.71 A, and placing the electrons symmetrically on the molecular axis, separated by 0.48 A, the 
experimental value of the H2 quadrupole moment is reproduced. From this model, we find ~ r 2 and $ ~ r 4 , for 
small changes in £ about r e . 

Regarding the orientational dependence, U' and U" are functions of the z-displacement of molecular center-of-mass, 
and of the polar (9) and azimuthal (</>) angles. U depends on z c , so the translational and rotational parts of U'(z c \ 9, <f>), 
U"(z c ; 9, cf>) do not decouple exactly. To separate the variables, we treat the interaction of translational and rotational 
modes in a mean-field manner [ p9[ . The effective interaction was computed by averaging the derivatives of Eqn.|5| 
over the center-of-mass vibrations of the molecule. Treating these as modes of a simple harmonic oscillator, the rms 
deviations in the z-position are found to be Sz rms — 0.298 A on a single graphene surface, and Sz rms — 0.225 A for 
H2 in an IC. Next, to compute < U' > T and < U" > T , variations were performed with respect to r in the adsorption 
potential for different r. Finally, we averaged over all orientations (#,</>) of the molecular axis. 

Using the values a=-1.6 and _B e /?ia; e =0.0138 pjj, we find that for the case of H2 adsorbed on graphene, our 
procedure predicts an upshift of Acji^o = 1.4 cm -1 (< iaU' >= —126.2 cm -1 , < U" >— —28.3 cm -1 ). In contrast, for 
H2 adsorbed in an IC, we predict a downshift of Awi^o=-2.9 cm -1 (< 3all' >= 168.0 cm -1 , < U" >= —42.2 cm -1 ). 
We expect similar results for D2: the same potential parameters and quadrupole moment may be assumed, though a 
and B e /hu} e are smaller |l4|], and should lead to smaller shifts. In the case of HD, the center of mass and charge do 
not coincide, leading to a larger Lennard- Jones contribution to the potential and a larger upshift. 
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In summary, we have observed several lines within the rotational and vibrational Raman spectra of H2, HD, and D 2 
adsorbed on SWNT; observed frequency shifts are small and consistent with physisorption. Comparative Q-branch 
data for H 2 on Cgo reveal substructure in the Q-branch, which appears to arise from H 2 adsorbed in different sites. An 
interaction potential was developed to estimate the frequency shifts of the adsorbate vibrational modes; our results 
support the multiple-site interpretation of our data. We wish to thank Vin Crespi, Roger Herman, John Lewis, and 
Dragan Stojkovic for their generous contributions to this work. This work was funded by the Department of Energy 
and the Army Research Office. PCE and KAW were funded, in part, by the Office of Naval Research. 



[1] M.S. Dresselhaus, K.A. Williams, and RC. Eklund, MRS Bulletin 24(11), (1999), 45-50. 
[2] K.A. Williams, and P.C. Eklund, Chem. Phys. Lett. 320(2000), 352-358. 

[3] M.M. Calbi, S.M. Gatica, M.J. Bo jan, G. Stan and MW. Cole, submitted to Rev. Mod. Phys.; 

http:/ /xxx. lanl.gov/abs/cond- mat/0103607 . 
[4] M.W. Cole, V.H. Crespi, G. Stan, J.M. Hartman, S. Moroni, and M. Boninsegni, Phys. Rev. Lett. 84, 3887 (2000). 
[5] A.C. Dillon, and M.J. Heben, Applied Physics A 72(2001), 133-142. 

[6] M.J. Heben, A.C. Dillon, T. Gennett, J.L. Alleman, P.A. Parilla, K.M. Jones, and G.L. Hornyak, Proc. Mat. Res. Soc., in 
press. 

[7] C.K.W. Adu, G.U. Sumanasekera, B.K. Pradhan, H.E. Romero, and P.C. Eklund, Chem. Phys. Lett. 337 31-35 (2001). 
[8] W.C. Conner, in Physical Adsorption: Experiment, Theory, and Applications, F. Fraissard (ed.), NATO ASI Series C 491. 
[9] AD. Buckingham, Proc. Roy. Soc. A 248, 169 (1959). 
[10] A. Kuznetsova, D.B. Macwhinney, V. Naumenko, J.T. Yates Jr., J. Liu, and R.E. Smalley, Chem. Phys. Lett. 321(2000), 
292-296. 

[11] S.A. Fitzgerald, T. Yildirim, L.J. Santodonato, D.A. Neumann, J.R.D. Copley, and J.J. Rush, Phys. Rev. B 60(9) 6439- 
6451 (1999). 

[12] J. Humlfcek, Journal of Quantitative Spectroscopy and Radiative Transfer, 27(4) 437-444 (1982). 
[13] AD. May, G. Varghese, J.C. Stryland, and H.L. Welsh, Can. J. Phys. 42(1964), 1058-1069. 

[14] G. Herzberg, Molecular Spectra and Molecular Structure: I. Spectra of Diatomic Molecules, (Van Nostrand Reinhold, New 
York, 1950). 

[15] K.P. Huber, G. Herzberg, Molecular Spectra and Molecular Structure: IV. Constants of Diatomic Molecules, (Van Nostrand 

Reinhold, New York, 1979). 
[16] G. Norlen, Physica Scripta 8 249-268 (1973). 

[17] K. Burns, K.B. Adams, and J. Longwell, J. Opt. Soc. Am. 40(6) 339-344 (1950). 

[18] CM. Brown, T. Yildirim, D.A. Neumann, M.J. Heben, T. Gennett, A.C. Dillon, J.L. Alleman, and J.E. Fischer, Chem. 

Phys. Lett. 329 (2000) 311-316. 
[19] L.W. Bruch, Surf. Sci. 125, 194 (1983). 
[20] H.R. Philipp, Phys. Rev. B 16, 2896 (1977). 

[21] W.B.J.M. Janssen, T.H.M. van den Berg, and A. van der Avoird, Phys. Rev. B 43, 5329 (1991). 

[22] F.Y. Hansen, L.W. Bruch, and S.E. Roosevelt, Phys. Rev. B 45, 11238 (1992). 

[23] H.Y. Kim, and M.W. Cole, Phys. Rev. B 35, 3990 (1987). 

[24] S.C. Wang, L. Senbetu, and C.W. Woo, J. Low Temp. Phys. 41, 611 (1980). 

[25] G. Maroulis and D.M. Bishop, Chem. Phys. Lett. 128, 462 (1986). 

[26] G. Stan and M.W. Cole, Surf. Set. 395, 280 (1998). 

[27] These calculations were performed using the PARATEC package developed at U.C. Berkeley. 

[28] R.M. Herman and S. Short, J. Chem. Phys. 48, 1266 (1968). 

[29] AD. Novaco and J. P. Wroblewski, Phys. Rev. B 39, 11364 (1989). 



4 



Adsorbate 


Transition 




Vads 


A 

L\V 


species 




(cm J 


(cm J 


(cm ) 


H 2 


Q(o) 


4160.5 


4162.3 


+1.8 








4161.1 


+U.0 








4ioy.o 


-l.z 




Q(i) 


4154.4 


4156.4 


+2.0 








4100.Z 


+U.0 








/II CO K 


-u.y 




b(0) 


9 C /I A 

O04.4 


OK/I Q 

oo4.o 


+U.4 




b(l) 


00/. U 


CQ7 r 

OO / .0 


+U.0 


HD 


Q(o) 


3629.8 


3632.5 


+2.7 




Q(l) 


OOZO.O 


OflOO Q 
ODZO.O 


+U.D 




rvo\ 
Q(2) 


ODzU.D 


ocon n 

oozu.y 


+U.4 




b(0) 


Oft? 1 

ZD / .1 


ZD / . / 


+U.0 




b(l) 


44o.l 


/I /I /I 1 
444.1 


1 A 
+ 1.0 


D 2 


Q(0) 


2993.5 


2993.8 


+0.3 




Q(i) 


2991.4 


2991.9 


+0.5 




Q(2) 


2987.2 


2987.5 


+0.3 




S(0) 


179.1 


178.5 


-0.6 




S(l) 


297.5 


297.8 


+0.3 



•5 




350 400 450 

Raman Shift (cnv 1 ) 

FIG. 1. S-branch data for Eh, HD, and D2 on SWNT at 8 atm, 5 atm, and 8 atm, respectively. Asterisks indicate unrelated 
lines due to external calibration sources. 
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FIG. 2. Q-branch data for H2, HD, and D2 on SWNT at 8 atm, 5 atm, and 8 atm, respectively. For comparisons of these 
peak positions with the frequencies of the free molecules, see Table 1. 
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